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Abstract Long range Rydberg blockade interactions 
have the potential for efficient implementation of quan- 
tum gates between multiple atoms. Here we present and 
analyze a protocol for implementation of a fc-atom con- 
trolled NOT (CfcNOT) neutral atom gate. This gate 
can be implemented using sequential or simultaneous 
addressing of the control atoms which requires only 
2k + 3 or 5 Rydberg it pulses respectively. A detailed 
error analysis relevant for implementations based on al- 
kali atom Rydberg states is provided which shows that 
gate errors less than 10% are possible for k = 35. 

Keywords Quantum computing • Rydberg atoms 



1 Introduction 

Quantum gates and entanglement of pairs of neutral Rb 
atoms have been demonstrated in the last few years us- 
ing Rydberg blockade interactions [TJUHllIS] . ^ e ^ ea 
of using transient excitation of Rydberg states to medi- 
ate strong interactions was introduced more than 10 
years ago [5] and has attracted a strong interest due 
to the potential for implementation of a wide variety 
of quantum computing and quantum communication 
protocols [7]. 

The fundamental Rydberg blockade two-qubit gate 
relies on the fact that excitation of a single atom to a 
Rydberg level can block the subsequent excitation of a 
second atom when the dipolar interaction between two 
Rydberg atoms is much stronger than the light-atom 
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coupling driving the excitation. If the Rydberg excita- 
tion lasers are resonant for transitions between one of 
the qubit levels and the Rydberg level, the excitation 
of the first atom is dependent on the qubit state, and in 
this way the excitation of the second atom also becomes 
conditional on the state of the first atom. Conditional 
evolution of the second atom is sufficient for creating 
two-atom entanglement and, together with single qubit 
operations, provides a library of gates that can be used 
for universal quantum computation. 

Although any desired unitary evolution of a multi- 
bit register can be implemented by combining sequences 
of the one- and two-qubit universal gates [8] such a pro- 
cedure is often inefficient. For example the Toffoli gate 
or C2NOT gate where the logical AND of two control 
atoms determines the transformation of the target atom 
state appears frequently in quantum algorithms. Unfor- 
tunately a single Toffoli gate requires at least six CNOT 
gates for its implemcntation[9] which increases the cir- 
cuit complexity. 

Rydberg blockade interactions have the very use- 
ful property that excitation of a single control atom 
can simultaneously block the excitation of many tar- 
get atoms, provided they all are close enough, that is 
within a so-called "blockade radius" [TD] , from the con- 
trol atom. This radius can be on the order of tens of 
microns which could reach many qubits in both 2D and 
3D optical dipole trap arrays. The blockade interaction 
has been utilized to devise a number of protocols for 
multi-bit encoding schemes [11] and quantum gates [l2l 
[T3lfT4lfT5] that provide direct routes to entanglement of 
several qubits at a time. 

In this contribution we show that the k control atom 
(CfeNOT) generalization of the Toffoli gate can be im- 
plemented efficiently using Rydberg blockade interac- 
tions between multiple atoms. Our most efficient im- 
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plcmcntation requires only 5 Rydberg ir pulses, inde- 
pendent of k, which can be contrasted with the best 
known circuit implementation of a CfcNOT which re- 
quires 32k — 120 elementary 1- and 2-qubit gates for 
k > 5 [16]. It has further been demonstrated that a 
C fe NOT requires at least 2k + 2 CNOT gates [I7| which 
shows that the 6 CNOT construction of the Toffoli gate 
in [5] is optimal. We have explicitly shown that our 
CfcNOT protocol can be used for efficient implementa- 
tion of Grover's quantum search [T8]. and we anticipate 
that it will also be useful for other quantum algorithms. 

In Sec. [2] we describe two different implementations 
of a Rydberg CfcNOT gate using sequential and simul- 
taneous addressing of the control atoms, and discuss the 
major sources of error during gate operation. In Sec. [3] 
we analyze the sequential addressing approach in de- 
tail and give analytical and numerical estimates for the 
gate errors using parameters relevant for implementa- 
tion with cold alkali atoms. The analysis is repeated 
for the simultaneous addressing approach in Sec. |H We 
conclude in Sec. [5] with a discussion of the results ob- 
tained. 



2 CfeNOT pulse sequences 

The CfeNOT gate can be implemented using Rydberg 
blockade in two ways as shown in Figure [1] If we are 
limited to using individual addressing and a single Ry- 
dberg state, then the gate is implemented by applying 
a series of Rydberg tt pulses to the control qubits in 
order, i.e. from 1 to k, between the |0) state and an ap- 
propriate Rydberg state |r). Then a sequence of three 
tt pulses to this Rydberg level is applied to the target 
qubit that swaps the target state as long as no Rydberg 
blockade is present, and finally the control qubits are 
returned to their original state with a phase shifted tt 
pulse to each control qubit in reverse order, i.e. k — > 1 
as shown in Figure QJi. The total number of Rydberg tt 
pulses is thus 2k + 3. This 2k + 3 pulse sequence was 
originally presented in [7] Sec. III.C.3. 

If we have the ability to excite to two different Ry- 
dberg levels and can simultaneously address the con- 
trol qubits then another possibility emerges. Here all 
of the control qubits are excited to a mutually non- 
interacting Rydberg level \s) followed by swapping the 
state of the target qubit via a Rydberg level |r) that is 
blockaded by any control qubit which happens to be in 
the |s) state. Finally the control qubits are returned to 
their initial state again using a phase shifted tt pulse. 
This implementation requires only five tt pulses inde- 
pendent of k. As we will see this is advantageous since 
the amount of time spent in the Rydberg level is re- 
duced at large k. However, it is difficult to find pairs 
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Fig. 1 CfcNOT pulse sequences, a) Sequentially addressed 
CfcNOT gate where the numbers label the order of the pulses 
applied to the control bits. For this sequence the first control 
qubit in the |0) state will block all the remaining pulses until 
it gets de-excited. This blockade prevents the target state swap. 
Therefore the target state is swapped only when all the control 
qubits are in the |1) state, b) Simultaneously addressed version. 
Here all the control bits are excited to a Rydberg state |s) that 
does not interact strongly with the same Rydberg state allowing 
any control qubit in the |0) state to be excited. Next the target 
state is swapped via a Rydberg state \r) which has a strong inter- 
action with the \s) state. Finally the control qubits are returned 
to the ground state. Therefore if any control qubit is in |0) the 
target state swap pulses are blockaded. 



of interacting and non-interacting states with sufficient 
asymmetry (B rs 3> B ss ) to obtain very high gate fideli- 
ties. For the Toffoli gate (k = 2) an alternative scheme 
using three different Rydberg levels is also possible, as 
was proposed in [T5] . 

For either approach we could replace the three pulses 
on the target atom with a single 2tt pulse to implement 
a CfcZ gate which would require a total of (2k + 2) 
or four tt pulses, and then obtain a C&NOT by adding 
Hadamard rotations to the target atom before and after 
the gate. Both approaches, using CZ or the controlled 
swap approach of Fig. [TJ were demonstrated for two- 
atom gates in [3]. 



3 Error estimates for sequential addressing 

In this section we develop detailed error estimates for 
the sequential addressing version of CfcNOT. When such 
a gate is used in a quantum computer the errors will 
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vary depending on the state which the gate operates on. 
Some input states will lead to higher errors than oth- 
ers. Since the number of states is exponentially large an 
analysis which calculates the error separately for each 
possible state quickly becomes unwieldy We therefore 
follow the simpler procedure of averaging over all possi- 
ble input states without regard to the spatial location of 
the information in the qubit array. The resulting error 
estimates should therefore be considered average errors: 
specific input states may have substantially higher er- 
rors. 

Furthermore we consider only the intrinsic errors 
due to the physical implementation of the gate. Addi- 
tional technical noise sources such as laser pulse fluctu- 
ations, finite temperature Doppler and atomic spatial 
position effects, and magnetic fluctuations are not con- 
sidered. These additional noise sources may be signifi- 
cant, particularly when using the gate to create entan- 
gled states [23, but their analysis lies outside the scope 
of this contribution. 

3.1 Intrinsic errors 

There are several intrinsic error sources for implemen- 
tation of the CfcNOT gates. These errors arise from the 
finite lifetime of the Rydberg states and undesired exci- 
tations due to insufficient blockade shift or off-resonant 
excitation from the other qubit level. To derive the in- 
trinsic errors we follow the methodology of [21117] and 
assume that the errors from each source are small and 
can be considered independently. Then each contribu- 
tion can be added together to obtain a total error. Fi- 
nally, to obtain the error for the gate, we average the 
errors over all possible computational basis input states. 

First we consider the individual addressing case shown 
in Fig. [IJi). We will calculate in turn the following er- 
rors: a) spontaneous emission from Rydberg states for 
control (c) and target (t) atoms, and b) pulse rotation 
errors from finite blockade shift or off-resonant excita- 
tion of the other qubit state. For error a) each control 
qubit undergoes a ir pulse between the |0) and Ryd- 
berg states, a wait period, and a final tt pulse back 
to the ground state. For each tt pulse the spontaneous 
emission error is given by ^7 where fl is the Rabi fre- 
quency between the ground and Rydberg state and r is 
the Rydberg state lifetime. During the waiting period 
the Rydberg state can decay giving an error of - where 
t is the time spent in the Rydberg state. In addition 
higher order terms occur when the qubit is excited by 
an off resonant pulse either because of finite blockade 
shifts or because the qubit was in the |1) state which is 
not connected to the Rydberg levels. The pulse scheme 
gives an error of ^7 + "^"^ + 5757 wnen the qubit 



is in the |0) state. In this formula n wa ;t is the number 
of 7r pulses that happen between the two pulses that 
target the control qubit i which can be found to be 
n wa it = 3 + 2 (fc — i) where the 3 comes from the 3 tt 
pulses applied to the target qubit and the other term 
accounts for the pulses applied to the i + 1 to k control 
qubits. Next we must count the number of input basis 
states that have this error and average this over the 
computational basis states. This error can only occur 
when the qubit i is in |0) and no earlier control qubits 
were in this state to cause Rydberg blockade. This num- 
ber of states is found to be 2 k+1 ~ l for control bit i and 
since there are 2 k+1 possible computational states we 
find that by averaging over the input states we get a 
factor of 1/2*. Putting this together we find that the 
leading order error term due to spontaneous emission 
for the control qubits is 

i=i 

There is an additional spontaneous emission error 
from off resonant excitation due to blockade leakage. 
For this case we do not have full excitation to the Ry- 
dberg state and the amount of population is given by 
(2 2 /2B 2 where B is the blockade shift of the Rydberg 
state. For this we assume that the blockade shift for 
all qubits can be given by some common minimum 
blockade shift and thus we will conservatively calcu- 
late the maximum possible error. The amount of error 
is the same as the resonantly excited case just multi- 
plied by the excitation probability of I2 2 /2B 2 ; however, 
the number of states seeing this error is different. Ac- 
counting for this difference we find 

_ nf2 ^ l + 3 + 2(fc-i) ^ ofe _. 

Use,c,2 - 2B2t 2^ 2 k+l 2^ Z 

i=2 j=l 



The spontaneous emission errors for the target qubit 
are a little easier to find as counting the number of 
states for each error is easier. The target sees a pulse 
scheme with three Rydberg tt pulses. First a pulse from 
|0) to Rydberg, then from |1), and finally from |0) again. 
We assume that the Rabi frequency for both ground 
states is the same. For the lowest order term we find 

F -—- 
and the blockade leakage term is found to be 

5ttJ? 2 fe - 1 



1 



The other type of errors b) can be termed rotation 
errors. These are due to off-resonant excitations which 
cause residual undcsircd rotations between the ground 
and Rydberg states. The leading term is due to block- 
ade leakage and is proportional to f2 2 /2B 2 . Since the 
blockade shift is generally small compared to the qubit 
frequency for implementations involving alkali hyper- 
fine states the term proportional to f2 2 /2ujf is smaller, 
where oj\o is the ground state splitting. Performing an 
analysis similar to that for the spontaneous emission 
errors we find 
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It is important to note that we have accounted for 
changes in the effective ground state splitting due to 
the blockade shift. In general the blockade shift for this 
off resonant state will be different from the blockade 
shift for the targeted ground state due to other nearby 
Rydberg levels but for simplicity we have assumed that 
it is still given by the same minimum blockade shift. 

Combining the eight errors listed above we get a 
total C fc NOT error of 
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Several comments can be made about the total error. 
The Rydberg state spontaneous emission error propor- 
tional to I/Qt enters with weight k due to the time 



needed to attempt excitation of all k control atoms, 
plus a term ~ l/2 fe from the target atom which is the 
small probability that the target atom is not blockaded. 
Thus the target atom error spontaneous emission error 
is negligible. We see that all terms except the first one 
scale linearly or slower with k. The error contribution 
which is quadratic in k comes from (TTJ) for the sponta- 
neous emission of the control bits. There are k control 
bits, and the amount of time an average bit spends 
in a Rydberg state is proportional to fc, hence the er- 
ror scales as k 2 . The prefactor multiplying this term is 
smaller by a factor of l/(J?r) <C 1 than the second line 
in Eq. ([2]) which implies that for moderate values of k 
the error grows linearly in k, after which it will switch 
over to quadratic in k for large k. 

The gate error is minimized by choosing the opti- 
mum value of J? with respect to the available block- 
ade shift B and Rydberg lifetime r. The relative scal- 
ings for implementation of Rydberg blockade gates are 
wio ^>B^>l2^>l/r. Dropping the terms ~ l/w 2 and 
working in the limit of Bt > fc we find the optimum 
Rabi frequency to be 



opt 
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which is the same scaling as for the CNOT gate[3T|. Nu- 
merical checks show that the optimum frequency found 
numerically agrees with J7 op t to within about 10% pro- 
vided Br > lOfc. 

Using i? opt we obtain an approximate error estimate 

of 
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r 4/3 



opt 
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(4) 



The (Br) 2 / 3 scaling is the same as that found for the 
CNOT gate[21]. For (Br) 2 / 3 > k the C fc NOT gate error 
grows linearly in fc. The less favorable quadratic in fc 
error scaling takes over when fc becomes comparable to 
(Br) 2 / 3 . 

Equations (|2l4j) are based on the assumption that 
all pairs of atoms have the same blockade interaction 
B. In a real implementation in an optical trap array or 
lattice the pairwise separations will range from some 
minimum to maximum value, and it is therefore nec- 
essary to take the spatial arrangement of atoms into 
account. We provide this analysis in Sec. 13.21 below. 

Nevertheless, if we simply assume that all atom pairs 
have the same blockade shift as the weakest interact- 
ing pair we will obtain some insight into the order of 
magnitude of the gate errors. As an example of this 
type of estimate let us consider Cs ns Rydberg states 
in a lattice where the largest two-atom separation is 
R = 20 fj,m. The C fc NOT errors for n = 100, 125, 150 
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Fig. 2 CfeNOT error assuming the same interaction strength 
for all atom pairs in the sequential addressing scheme for Cs 
ns states with n = 100, 125, 150. Parameters used were r = 
(330,540,820) fia, B/2tt = (0.69,9.0,52.) MHz, and u) 10 /2n = 
9200 MHz. The blockade shift values were calculated for a fixed 
two-atom separation of R = 20 ^m. With the above parameters 
we find ft opt /2-7r = (0.11,0.53, 1.5) MHz. 

are shown in Fig. [2] The displayed curves were calcu- 
lated using the approximate result Q as well as from 
([2]) with Q — > f2 opt . On the scale of the figure there is no 
discernible difference between the two error estimates. 
We see that the error decreases rapidly with principal 
quantum number and is about 0.06 at k = 50 for the 
150s state. The error grows quadratically at much larger 
values of k than shown in the figure. We emphasize 
that the errors shown in the figure are only meant to 
be thought of qualitatively since they assume all atom 
pairs have the blockade shift of the pair with the largest 
separation. In order to find a more accurate error esti- 
mate we must take into account the physical placement 
of the qubits as is done in the next section. 

3.2 Lattice averaging 

We will assume that the atoms are placed on a two- 
dimensional square lattice with period d, as shown in 
Fig. [21 While a three dimensional system would give 
a larger number of atoms within a blockade sphere, it 
appears difficult to achieve individual control of each 
atom in such a geometry. We will therefore concentrate 
on a two-dimensional implementation. 

There are several factors that limit the minimum 
value of d that may be used. In order to avoid collisions 
between a Rydberg excited electron and a neighbor- 
ing ground state atom we require d >^ a^n 2 with an 
the Bohr radius. A precise evaluation of the maximum 
number of atoms that can be placed within a blockade 
sphere while meeting this requirement can be found in 
[2"2"] where it was shown that in 2D k' ~ n 2 l z at con- 
stant gate error (we use k' = k + 1 for the number of 
atoms) . 
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Fig. 3 Two-dimensional qubit lattice with period d. The central 
atom serves as the target qubit. 

This scaling needs to be modified by the additional 
limitation that B(d), the blockade shift between atoms 
separated by d, should satisfy B(d) < \U(n) — U(n — 
l)]/2 ~ n~ 3 , where U (n) is the energy of Rydberg level 
n. If B(d) exceeds this amount then the Forster reso- 
nance between pairs of atoms with n differing by ±1 
may lead to an interaction induced resonance, which 
would break the blockade effect. Setting B(d min ) = [U(ri)— 
U(n — l)]/2 with resonant Forster scaling B(d) ~ n 4 /d 3 
gives d m ; n ~ n 7 / 3 . This is slightly more restrictive than 
the collisional limit in the preceding paragraph. When 
the largest two-atom separations are in the van der 
Waals regime the blockade scaling is B(i?) ~ n 12 /R 6 for 
the heavy alkalis [33] which leads to an R limit at con- 
stant gate error of i? max ~ rJ I 3 . The number of atoms 
scales as k' ~ (i?max/d m in) 2 ~ constant. We conclude 
that, asymptotically for large n, the number k' is con- 
stant and there is no advantage to working with very 
large values of n. Indeed, the numerical results given 
below show that the optimum performance is obtained 
for n ~ 75. 

The reason for this preference for moderate values of 
n is somewhat subtle and can be understood as follows. 
The Forster interaction B(n,n') is different between an 
atom pair (n, n) and an atom pair (n, n' ^ n). Calcula- 
tions show that for ns states with n > 100 the interac- 
tion varies smoothly with n'—n and B(n, n) > B(n, n' =/= 
n) for \n' — n\ small. However, the difference in interac- 
tion strength for n' = n — 1 is only about 20% so the 
limit on d m ; n indicated above could only be slightly re- 
laxed. This situation is qualitatively different at smaller 
n since the Cs ns Forster defects increase rapidly for 
n <^ 80.|23j. It turns out that the limit on d m i n is then 
no longer set by the need to avoid interaction induced 
resonances, but by limits on optical beam focusing, and 
the need to avoid crosstalk between neighboring sites. 
This limit is not fundamental and can be improved on 
by using beams with specific spatial profiles. For defi- 
nitcness we will assume a lower limit of d m ; n = 1 /im. 
With a fixed d m i n there are diminishing returns in going 
to lower n since the interaction gets too weak. We find 
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Fig. 4 Blockade shift for several Cs ns levels. The blockade is 
calculated as described in I23I . 
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Table 1 Sequential addressing CfcNOT gate errors from Eq. J2j 
averaged over a square lattice for several different Rydberg Cs 
ns levels. For n = 75 the Rabi frequency was set to f2/2n = 
(45, 29, 20, 15, 12) MHz for k = (3, 8, 15, 24, 35). 



numerically that the optimum tradeoff is in the vicinity 
of 7i = 75. 

We present representative numerical results for the 
gate error for Cs 50s — 150s states in Table [JJ The 
errors are found by averaging Eq. (0) over all values of 
Rij = |Ri — Rj| that occur in a square lattice, using 
the B(i?) curves from Fig. [4] This gives an error which 
depends on Q which we then set to minimize the error. 
In this way we capture both the resonant dipolc-dipolc 
and van der Waals behavior. For n = 50, 75 we set 
d = 1. /im. For n > 100 we improve on the resonant 
interaction limit described above by reducing d such 
that B(d) = 1.5 x [U (n) — U(n — 1)]. In a square lattice 
the next separation is R = \[2d which gives a blockade 
shift of B(d)/2 3 / 2 ~ 0.5 x [U(n) - U(n - 1)]. Further 
optimization with B(d) several times larger than the 
gap may be possible but this requires additional fine 
tuning of the architectural parameters which is a task 
we leave for future work|24j. 

We see that for k < 8 the error is approximately 
constant with n whereas for k > 15 the 75s level gives 
the best results. We obtain an intrinsic gate error less 
than 0.1 for k = 35. The variation of the error with 
Q is shown in Fig. [5] for k = 24. At lower n the Ryd- 
berg lifetime is shorter and the optimum Rabi frequency 
is about 8 MHz, decreasing down to about 1 MHz at 
n = 150. The numerically found optima are about 20% 
lower than predicted by Eq. ((H) and they occur at about 
twice smaller Rabi frequencies than predicted by Eq. 



0.10 
0.08 



CD 

TO 0.04 
0.02 















125\ ' 




_ ~~ + 


100 


n=75 









2 4 6 8 10 12 14 

Q/2n (MHz) 

Fig. 5 Dependence of the gate error on J? for k = 24 and several 
Cs ns levels. The crosses near each curve mark the values of i? pt 
and Eopt found from lattice averaging of Eqs. II 3 1 1 1 > . 

• The difference between the analytical estimates and 
numerics can be attributed to lattice averaging of the 
interaction in the numerical solution. The error minima 
are sufficiently broad that no unreasonable fine tuning 
of f2 is needed. 

We emphasize that the estimated gate errors repre- 
sent averages over states of the k + 1 bits in the compu- 
tational basis. Specific input states will have gate errors 
that are larger than the average by some factor. In order 
to assess the strength of our results it is desirable to put 
some bound on the error variations. We first point out 
that our error estimates are conservative since we es- 
sentially compute the probability that "something goes 
wrong", and assign zero fidelity to the resulting state 
(even if it may actually be an algorithmically useful 
outcome despite the error). 

To be more quantitative we can estimate the er- 
ror for what is plausibly a worst case input state. Let 
us suppose a single control bit is Rydberg excited, but 
that it appears first in the Rydberg excitation sequence, 
and thereby spends the longest possible time in the Ry- 
dberg state, and that it has the largest possible spatial 
separation from the target bit. To estimate the error 
for such a case we use the parameters from Table [TJ 
n = 75 and use Eq. (26) from [7] for the gate error. 
We calculate the blockade shift B from Fig 2] for the 
largest two atom separation that occurs for k = 35 and 
artificially decrease the Rydberg lifetime r by putting 
r — > t/35 to account for the longer time spent in the 
Rydberg state. Doing all this gives a C2NOT error of 
0.02 which is still less than the average error in Table 
[TJ This may indicate that our estimate of the average 
error is, indeed, a conservative one. An identification of 
what are the worst possible input states, and what is 
the best choice of parameters to reduce the dominant 
errors remains an open question for future work. 

If we were to implement the C35NOT gate with 
elementary gates we would need at least 72 CNOT 
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gates [T7| . It is known from random benchmarking studies (2"5] 
that the total error grows approximately linearly with 
the number of gates, so an equivalent error for the 
C35NOT gate would require CNOT gates each with an 
error of 0.001 which is well beyond the current state of 
the art for gate error which has been achieved in ion 
traps PB]. 

As can be seen from Table Q] the gate errors grow 
faster than linearly in k which differs from the result of 
Eq. Q . There are two reasons for this discrepancy. First 
the ratio of /c/(Bt) 2 / 3 is not very small for the largest 
k values considered so there is some quadratic scaling. 
More significantly, since there is a very large variation 
of B across the array it is not possible to simultaneously 
optimize all errors with a single fl value, which violates 
the assumption leading to fl3J). 



4 Error estimates for simultaneous addressing 

The simultaneous addressing scheme of Fig QJ>) can 
be analyzed along similar lines. Parameters associated 
with the control and target qubits are labeled with sub- 
scripts c and t, respectively. We identify the follow- 
ing five leading error sources due to spontaneous emis- 
sion and pulse rotation errors on the control and target 
qubits: 

-k + ^^k, (5) 
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The binomial sums that appear in several terms ac- 
count for averaging over multiple Rydberg excited con- 
trol atoms. The blockade shift between control and tar- 
get atoms is B ct and D cc is the direct interaction be- 
tween two control atoms which we calculate as the dipolc- 
dipole interaction strength between the laser excited 
states. Parameters will be chosen to ensure that 
|D cc /^cc| < 1 so there is no blockade between control 
atoms. We assume the blockade and interaction shifts 



N 

1 



Q. 

m 



J 000 




Fig. 6 The blockade B c t and interaction D cc strengths for con- 
trol atom in \s) = |60p3/2> m = 1/2) and target atom in 
\r) = \60s 1 /2,'fn = 1/2} states. 



are additive in the number of Rydberg excited control 
atoms, which can be as large as k. For k large the is- 
sue of interaction induced Rydberg resonances which 
were discussed in the previous section may again be 
problematic, as well as nonlinear behavior of multiple 
interacting Rydberg atoms [27] • While detailed study of 
such effects is outside the scope of this work, we will see 
that the simultaneous addressing version of the gate can 
be operated in a regime where B ct (<i), D cc (<i) are small 
compared to the Rydberg gap U(n) — U(n — 1) which 
adds confidence to the validity of our present results. 
Note that fl c is only constrained by the off-resonant 
excitation limit of i? 2 /u; 2 Q so we expect that the mini- 
mum error will be found at large values of fi c . 

To proceed we add Eqs. ((5j - ^ to arrive at a to- 
tal error which is averaged over the lattice. We then 
minimize the error with respect to both f2 c and fi t - 
In the averaging procedure D cc is only averaged over 
the Rij between control atom sites and B ct is only av- 
eraged over Ri j t where jt labels the target atom po- 
sition in the center of the array and i runs over the 
control atoms. The |s) and \r) Rydberg levels used for 
control and target qubits respectively should be cho- 
sen to maximize B ct and simultaneously minimize D cc . 
This type of asymmetric Rydberg interaction has been 
invoked previously for multi-bit gates in [T9l[14"] . We 
choose the control-target interaction to be dipole al- 
lowed so B ct ~ n 4 /R 3 and seek parameters for which 
the control-control interaction is in the van der Waals 
limit where D cc ~ n 12 /R 6 . The asymmetry thus scales 
as B ct /D cc ~ R 3 /n s which suggests using relatively 
small 71. Numerical search through different Rydberg 
levels and lattice periods yielded acceptable results for 
\s) = |60p 3/2 , m = 1/2) and |r) = |60s 1/2 , m = 1/2). In- 
teraction curves for these states are given in Fig. [6] Note 
that since we are now using p states the interaction is 
no longer isotropic and a quantization direction has to 
be defined. In order to avoid having to perform lattice 
averaging of the angular dependence we have assumed 
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d (fJ.ro) 


k = 3 


8 


15 


24 


35 


error (77 K) 


4.2 


.02 


.03 


.07 


.09 


.16 


error (300 K) 


4.0 


.03 


.04 


.10 


.14 


.23 


H c /2n (MHz) 


4.0 


180 


240 


300 


350 


390 


flt/2-K (MHz) 


4.0 


0.81 


1.6 


1.3 


1.8 


1.6 



Table 2 Simultaneous addressing C^NOT gate errors averaged 
over a square lattice. Rydberg states were \s) = |60p3/2i m = 
1/2), \r) = \60s 1/2 ,m = 1/2) with lifetimes taken to be[28l 
T 60p 3 / 2 = (335,148), T60s 1//2 = (171,97) (j,s at temperatures of 
(77K,300K). The Rabi frequencies listed are for the 300 K case. 



the quantization axis is perpendicular to the plane of 
the array so the dipole-dipole interaction is isotropic in 
the plane. This may not be the optimum orientation 
as far as the gate errors are concerned, so the results 
obtained need not represent the best possible solution. 

Results for the gate errors and corresponding Rabi 
frequencies are given in Tabled For the states used and 
an optimal lattice period of about 4.0 /jm, the Rydberg 
interaction strengths are a small fraction of the Ryd- 
berg level spacing so interaction induced resonances are 
not an issue for this approach. We see that the errors 
scale slightly slower than linearly with k, but they are 
about three times larger than for the sequential scheme 
at k = 35. The errors can be reduced by about 30% 
by going to liquid nitrogen temperatures (first line of 
the table) which increases the Rydberg state radiative 
lifetimes. 

Although the errors are larger it is worth noting 
that the number of pulses is about 15 times smaller 
compared to the sequential approach for k = 35. In ad- 
dition the gate time is T = Sir/ Q t + 2ir/f2 c = 1.1 fis for 
k = 35 compared to about 9 //s for the sequential gate. 
These differences imply better immunity to technical 
noise sources, the presence of which will influence the 
preferred approach in practice. 

5 Conclusion 

We have presented two implementations of a CfcNOT 
gate using Rydberg blockade interactions. Using realis- 
tic numbers we envision implementation of these gates 
with k as large as several tens of atoms and errors of 
only a few percent. The errors are optimized for k > 15 
at moderate values of n ~ 75 for sequential addressing 
and n ~ 60 for simultaneous addressing. The sequen- 
tial scheme has lower intrinsic errors, while simultane- 
ous addressing requires fewer laser pulses and is faster, 
which suggests a lower sensitivity to technical errors. 
These results may be contrasted with the situation for 
two-qubit gates where higher values of n are necessary 
for minimum errors at long range. The preference for 
smaller n for this multi-qubit gate stems from the need 



to pack as many bits as possible into the strongly inter- 
acting resonant dipole-dipole regime. This can be done 
more efficiently at smaller n values. 

As we have shown in [18] a slight modification of 
these CfcNOT gates can be used to efficiently imple- 
ment Grover search. The error analysis for the modified 
gates is slightly different due to the absence of the tar- 
get atom in the sequential addressing version. For com- 
pleteness we give the corresponding error expressions 
in the Appendix. The errors for the modified gates will 
be comparable to, but slightly less than, those in Table 
[TJ Although the CfcNOT gates cannot be scaled to ar- 
bitrarily large k with acceptable errors the sub-register 
architecture in [18] presents a path to scaling of the 
search algorithm to arbitrarily large search problems. 
The sub-register architecture could in principle also be 
used to scale our CfcNOT gate to larger values of k. 

The numerical analysis of errors has used Cs atom 
parameters. We expect similar results would be ob- 
tained for Rb, and possibly for other species being con- 
sidered for Rydberg gate experiments. The results pre- 
sented use realistic experimental parameters but could 
be further optimized along several directions. We have 
considered a controlled amplitude swap type of opera- 
tion which requires three Rydberg n pulses on the tar- 
get atom. This could be replaced by two singlc-qubit 
Hadamards on the target bracketing two Rydberg tt 
pulses, which give a CfcZ, to give a CfcNOT gate. The 
Rydberg state related errors on the target qubit would 
thereby be multiplied by a factor of 2/3. The Hadamard 
pulses could also be dispensed with, while keeping only 
two Rydberg pulses, using the bright /dark state cou- 
pling mechanism from|29j. In the simultaneous address- 
ing approach the gate fidelity is limited by the available 
asymmetry between D cc and B ct . One possible route to 
improving the performance would be to consider mi- 
crowave field dressing of the D cc Forster interaction. Al- 
though microwaves have mainly been viewed as a means 
of strengthening blockade interactions [3D], in principle 
it may also be possible to use them to increase the level 
shift of neighboring states, and thereby reduce the in- 
teraction strength. Furthermore, we anticipate that de- 
velopment of tailored pulse shapes or composite pulses 
will be a fruitful direction for future improvement of 
the gate errors in both approaches. 

In addition to the use of the modified C^NOT gate 
for quantum search we expect this gate will be useful 
for other quantum algorithms and for quantum error 
correction. In light of the high cost of implementing 
the CfcNOT using elementary 1- and 2-bit gates our 
approach may prove advantageous for building multi- 
qubit processors. Analysis of models for quantum com- 
putation predict error thresholds for fault tolerant oper- 
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ation at 1%[3T] or higher [32] • For large values of k our 
predicted errors approach 10%, and therefore exceed 
known fault tolerance thresholds. However, it should 
be remembered that the calculated thresholds assume 
1- and 2-bit gate operations, many of which would be 
needed to implement a single CfcNOT operation. It there- 
fore remains an open question as to the usefulness of 
many-bit gates with moderate error levels for fault tol- 
erant implementations of quantum algorithms. 
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6 Appendix 

Here we give the error estimates for the sequentially ad- 
dressed modified CfcNOT operations used in [IS]. The 
analysis follows that in Sec. [3] with the only difference 
being that we apply the 2k it pulses to the "control" 
atoms, but there is no target atom. In this way a con- 
ditional phase is picked up on the qubit register which 
is used to implement the Grover iterations. The error 
estimates are therefore slightly smaller and are given 

by 



E 



se.c.l — 



E 



se,c,2 



E r,c,l ~ TTi^r 



E r ,c,2 — 



7T 

4B 2 t 
2B 2 

a 1 



2fc-3 + 
k 2 - Ak 
k-2+- 



r 6 



'in 



1 - 



1 



fc-i 



O 2 



(wio ± sr 



1 



- 1 



The combined gate error is 
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In the limit of Br ^s> k ^> 1 the optimum Rabi fre- 
quency (2 op t and the minimum gate error are the same 
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as in Eqs. (|3I4[) . There is no difference since in the limit 
of k 1 the addition of the three target atom pulses 
in the CfcNOT protocol compared to a Grover iteration 
step has negligible impact on the error. 



